S  2 


NPSEC-93-005 


NAVAL  POSTGRADUATE  SCHOOL 

Monterey,  California 


New  Algorithms  for  the  Detection  and  Elimination 

of  Sine  Waves  and  Other  Narrow-Band  Signals  in  the 

Presence  of  Broadband  Signals  and  Noise 

by 

Michael  A.  Soderstrand 
University  of  California,  Davis 

K.  V.  Rangarao  and  H.  H.  Loomis 
Naval  Postgraduate  School 

October  1992 


Distribution  unlimited. 


FedDocs 

D   208.14/2 

NPS-EC-93-005 


Naval  Postgraduate  School 

Monterey,  California  93943-5000 


Rear  Admiral   R.W.    West,    Jr.  H.  Schull 

Superintendent  Provost 

This  work  reported  herein  was  partially  supported  by  the  Secretary  of  the  Air  Force. 

This  report  was  prepared  by  Michael  Soderstrand,  K.  V.  Rangarao,  and  H.  H.  Loomis, 
Jr. 

Submitted  by:  Reviewed  by: 


JJncjass  if  fed 


'-  ^"TSECdRi.vcLAss^cAno- 
Unclass  if  fed 


REPORT  DOCUMENTATION  PAGE 


NaJIT  V*  :ERf°^"^  organ^atio" 
Naval    Postgraduate    School 

Monterey,    CA  S3943-5000 


ea    NAME  Of  t-UIOiNo,  SPONSOR^; 
ORGANIZATION  ur^OR,l 


|  bo    OrFiCE   SYMBOl" 
(If  applicable) 

Code    EC 


D    RESTF       ive  Mark,\gs~ 

"-  ^auT,oN  B.,OMEPO- 

D'strfbutfon    unlimited. 

J       |,",Ula'0K!NG0^«...M,^^.^-r,>.V3£;J- 

Naval  Postgraduate  School 
Monterey,  CA  939^3-5000 


form  Approved 
OMB  Mo  0704-0  j 88 


8d    OFF  CE  s-V 
Of  applicable) 

ACBMB 


9P"       — -~~~~--     ^.•.-CA-oTaTaE?" 


pi     :  .  a  7" 

ELEME'v"   mo 


no 


ECT 


r  a  s  i 

NO 


NEW  ALGORITHMS    FOP    THE    OFTFrr.oM  J [ 

fm^y^^^^  "«»  *»  om*  NARR0W.BAND 

^■chael    A.    ^a.,t  ■ 

13a    ivpe  of  repor'T f       '     V-     Ranoaran,    Jnd    ^  H       ,  ^m?_ 


WORK    UN'T 

accession  no 


FlE.O 


_co: 

GROl 


^T    CODES 


SUB-GPO 


Subject 


ERVS  (tcnr/DLe  Cn  reverse    f" 

narrow-band    fnterfpr.n.  """'*  and  ,den,,t>  b>  t»°«  "<™~ 

-,_.       ..  ' " LC'  rerence    removal 

Four  different  classes  of  adaptive  skna   cZ  7       """.  " ■ 1 

^tu::t^onb  :ip:Lr  ctrsffffnd  (4)  y*£S3£  tt°;  :x  neomt  ffi,ters 

■in*  a  new  ^^S^Z^^^T^  T^*  WUh°Ut  d^adinS  Performance 
_^MBPSP.)  signals  with  and  without  background  noLe  ,nterf«ence  ,n  broadband  Bi-Phase| 

r>KTPio  ,t.^..  


IC  USERS 


ABSTRACT  SFf   iP  T  -   r    a^  • 

"nclassf f JPn 


'itions  are  obr.   /i  •< 

S/N  0102-LF-014-6603 


CLASS^f  :  .'.  •    i,     ,     -,.      ,  A 


•^^TTSISSIF.^T.0N  OF  THIS  PAGE_ 


DD  Form  1473.  JUN  86  (Reverse) 


NEW  ALGORITHMS  FOR  THF 
DETECTION  AND  ELIMINATION  OF  SINE   WAVES 

IN  THE  PRESENrpTHnp  T^**™   SI^S    ^ 

THE  PRESENCE  OF  BROADBAND  SIGNALS  AND  NOISE1 

Naval  Postgraduate  School  Technical  Report 


Electrical  p„„   "iChaeI  A"  s°°erstrand 

^^T^<c°rter  srience  e 

y  ot  Laiif . ,  Davis,  CA   9561 


Department 
6 


ElectriJ^  *a"garao  and  H-H-  Loomis 

NaTal  C°:fe7^~0  Department 
Naval  Postgraduate  School 

Monterey,  CA   93943-5100 


ABSTRACT 

Four  different  classes  of  srtjnH, 
"*ed  to  eliminate  narrow-band  fnterfer  ^"f    Cance1^  can  be 
"gnal:   (1)  cascaded  second-orner  1  T!-  r™  S  br°adband 
'"-line  notch  filters,  (   1  o   n  '  .    ter5'  <2)  nidh-order 
and  (4)  high-order  ba  dpa ss      ~c-ncl,    ""T   n°i5e  "Meiers, 
structure  based  on  second-nrZ  !  c°ncelei"s-  Of    the  four,  a 
cancelers  is  found  to  °  rforf bettef  tT  'JJ""  L^d  *S  ^".l 
The  adaptive  algorithm  for  these     tt      \ ^  °ther  structures. 

'-an  and  Martin  and  modified  by  PetranHa  Mf^"  Pr°P°Sed  b* 

■  recraglia,  Mitra  and  Szczupaf. 

The  Kwan  and  Martin  ctn,M, 
complexity  without  dearadina  «« f  ""  be  reduCed  in  h^ware 
algorithm  that  out-per »   9.  y  f Th""^*1"9  *  ne"  *daP"^ 
algorithms.   This  new  structure  °=  the  °ther  known  structures  or 
elimination  of  narrow-ban   inter £"" rt'CUlar1^  »""«(  to  the 
BnHt-Key  (BPSK)  signals  l^^",,^-  S3:" 


work  was  s 


upported  in  part  by  the  United  Stat 


es  Air  Force. 


0.0  Introduction 


Phase  Shan-Key  CBPSK)  moduS  Jn      , V  a  5eread-spectrum  Bi- 
narrow-band  interference   Th«=   9  '  ln  the  Presence  of 
less  energy  than  the  broaden sig^l  ^V"^' ^  »*   h-e 
concentrated  over  a  narrow  bandwidth  fs     be"USe  2t  ls 
signal.   In  the  specie  JrttZtJTrl   "T  °Ut  '""  br°*^nd 
5  tD  ruminate  only  the  narrow-band  i„   '"tested  in,  the  goal 
processing  will  separate  broad-hanH       ' BrBn"-      Later 
signal.   However,  this  later   r     n°1Se  irom   the  desired 
-less  all  narrow-band  int^^S'i"'^"  T   "^  ^^ 
broad-band  processing  is  begun        eliminated  before  the 

The  most  common  approach  tn  =i ■  • 
interference  is  to  make'use  of  a  casrart   ?   "T   na—°and 
Based  on  extensive  analysis  Zh      Cascade  oi   notch  filters  Cl-8] 
Martin  developed  a  par  c  ,larlv  ^T"4*1  WDr''  Kwan  *"" 
Edified  by  Petraglia,  Mitra  Inn  fKStrUCtUre  C53  whlch  «« 
good  success.   However,  the  hardware  r^  l?'B'   "nd  aDplied  «"h 
performance  with  more  than two  r    re1"»-ed  to  achieve  this 
becomes  too  complex  for  easy  LZ   thre\lnt^^ing   signals 
high  sampling  rate  is  required  '  0Z     <T'    particu^rly  when 
of  this  report  is  to  introduce"   *■    the  maJOr  contributions 
that  maintains  the  performance  of  the  ^  IT""""  -month. 
requires  substantially  lesshfrri!     °rl9>nal  algorithm  but 
1-d  interference  when  there'™!  J^nt^^'  ™" 

o.  --w^n^^^^  * 

describes  four  different  structures  f^  signals.   Sections  1 
cancellation.   Section  2  compares  the  f    PtlVB   lnt^rence 
basis  for  the  selection  of  the  second  n°T   T^*   and  9ives  the 
the  best  design.   Section  3  gi  ves  detaT^T  band"PaS£  canceler  as 
algorithm  for  the  second-order  band  nVJ        lnformation  on  the  new 
experimental  data  to  demonstrate  its  ad!  T"1^  a"d  ^nsive 
algorithms  including  the  afoorith™  h  aovantages  over  other 
gives  detailed  examples  of  the  n    f  "^  and  Martln-   Section  4 
cancelling  narrow-band  BPSr  sinn^  ^ithm's  Per<°™ance  in 
BPSK  signals  and  noise.   Con usionY^  the.presen«  °<   broad-band 
including  a  discussion  of  the  hard       ^  ln  Section  5 
—   Section  6  contains'tnTr^ ^Tt^r^™ 
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1.0  Adaptive  Notch  Filters 
Notch  mters  for  re.ov.no  ^t^^rSSJS*^^ 

Thus  IIR  pipelining  may  become  an  important  issue. 

1.1   Second-Order  Cascaded  Notch  Filters 

The  second-order  notch  filter  is  used  in  cascade  and  in-line 
With  the  siqnal  as  shown  in  Figure  1-la.   The  transfer  function 

for  the  notch  filter  is  given  by: 


2(2-k2-k12) 


1  - 


n-i.-  2-k 

i.  r-.  o 


HN<^  =  " TT    ;   2,   -i     »   ■   —2 


(1-1) 


If  k,  is 


1  -  (2-k2-k12)z_1  +  (l-k2)z 

For  arbitrary  values  of  k,  and  k2,  this  is  a  symmetric  notch 
filter  with  unity  qain  at  DC  and  the  Nyquist  frequency.   If 

e  tconstant,  then  the  3db  notch  width  is  -l"^-^- 
Thus  k,  may  be  adapted  to  remove  one  "arrow-band  s  gnal    A 
cascade  of  such  filters  can  be  used  to  remove  multiple  narrow 
band  signals. 

1.2  Second-Order  Cascaded  Signal  Cancel er 

The  cascaded  second-order  signal  canceler  approach  shown  in 
Fioure i-ib  has  the  advantage  that  the  desired  signal  does  not 
plssthrouoh  the  adaptive  niter.   Instead,  the  .band-pass  niter 
is  u=ed  to  detect  the  narrow-band  signal  which  is  then  ^ubtractec 

™  the  desired  signal.   A  constant  3db  bandwidth  notch  can  be 
achieved  by  selecting  a  band-pass  niter  with  the  transfer 
function: 


y.  (n) 


<*>   Cascade  o,  Second-Order  Notch  Filt 


ers 


»>      ^scade  o,   Second-Order  CanceHation  Filter, 


x(n) 


y(n) 


y(n) 


<c)      Higher-Order    Notch    Filter 


y(n) 


x  (n) 


(d)   Hi  qher-OrdPr  cln,  ■,  r, 

°er  Sl5nal  Cancelling  Filt 


y(n) 


er 


"■-.  1-1.   Four  Categories  of  Notch  Fl]t 


ers 


HEP(z)  = 


where:   D(z)  is  the  same 


(1-2) 
2D(z) 

denominator  as  HM(z)  in  equation  (1) 


The  signal-canceler  structure  is  also  ^^g^tSn. 

because  it  is  relatively  ^^"^X^t  to  the 
which  are  related  to  the  gradient  o^H„jjith  r^^  ^^ 

T^^ET^i"  oi ^his  'liter.   The  sensitivity  .unction 

Hs(z)=HEP<z)*H55(z>  where: 


2k i  z_i  d-3) 


H55(z)  = 


D(z) 
The  parameter  kl  .ay  then  be  adapted  by  the  formula:  I 

v  ,  v  (3) 

k±(n+D  =  ^i(n)  '  »e(n)s(n) 

Considerable  information  is  available  in_the  1"-^  « 

application. 

1.3  Higher-Order  In-Line  Notch  Filter 

r  r--    D  i_ir  hA=  the  advantage  that  it  is 
The  notch  filter "  of  F.gure  1  c  has  the    ^^^  as  , 

more  easily  pipelined.   S"7h.a-I  ««  t  Dlpeilne  and  has  the 
Hnear-phase  FIR  filter  which  is  easy  to   P^  ire 

advantaqe  of  linear  phase.   However,  an  FIR 
many  weights  to  obtain  good  performance. 

:r,rxi:  kiss  *«"» «'-  >• — •""""  "- 
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pipelining    the   FIR   -fiif 

^ccnd-order   IIR  fil\^   b£  {"•  -»«cult  than  pipelining  tne 

be  adapted  by  having  a  discrete  set   T        "'      ThiS  fUter  "uld 
could  be  switched   in   by  a  detectfnn  pre"sel^ted  notches  which 

««**,.   in   the  spectrl  of    the     nit    «CUlt,*hat   S'mpIy   looks     - 
approach  could  be  used  with   any  of   th      ]9?fK      filth°^h   this 
particularly  attractive  with  this  fnf    "J^5'    "    ls 
9-erate   geometric   pole/zero  patterns   fn   I "V   ^   """'V   to 

erns   ln  higher-order   filters. 

1.4        Higher-Order   Signal    Cancel er 

discussed  above.      However,    sLe  the^f'0,  ^   ^"^  notch 

the   signal,    the  desired   sinn*?  w  can"ler    is  not   in-line  with 

filter.      Also,    it        ^   -^"al    d„es  ™t   need   to  oass  throuoh   the 


1.5 


IIR  vs  FIR  Realizations 


notch  of  sectfon  .?aand  thense  ^r'^6""'  the  — d  order 
^ct  on  1.2  would  by  neces=ft   h    "^'^  band"Pass  cancel  er  of 
second-order  FIR  fiUer  cannot  T   '°  '*  HR  realizations   A 
band-width  in  order  to  ad^ff '^V  ^fi^tly   narrow 
Theoretically,  the  hiqher^orde rJlntj™*   lnte^-ang  siqnals. 
■4  could  be  either  IIR  Dr  pTo   **   nations  of  sections  1.3  and 
Nearly  suited  for  our JTU  tf'     ^ '  the  «cond-order  IIR  f« 
resonance  with  only  one  parget"  "E*  "  Pr°Vldes  a  sbarp     1S 
frequency,  to  be  adapted    Hioh!'  f   n°tch  °r   band-pass 
°reat  explication  without  s gnfficantl  ""  'lIUrS  "°Uld  ad° 

significantly  improving  performance. 

rrc  '""*"*  applications,  however   it  „   , 

FIR  rather  than  IIR  filters   fir',  ,/Sy  be  a°vantageous  to  us» 
»rder  to  realize  the  narrowband \tJrlV   ""  be  higher  °rder  „ 
primary  advantaqes  of  the  FIR  rl,  haracteristi".   The  two 
made  to  be  linear  phase  and  tha """"'."  that  the*  "n  be 
Discrete  Fourier  Transform^  such"  as".  F  ^  USin*  the 
"TT>  chip  or  a  Recursi    D    ''^"  p  aS„a  Fast  Fourier  Transform 

e  m-T  using  Residue  Number  Arithmetic 
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[9,   In  th.  instance  * F  -£!« ,  J  .3  an  N-U --  ^ nfDf  the 
adaptive  in-line  notch  filter  woul a  o      y        adaptive 
many  FIR  design  techniques   This  could  .«pl y  ^  ^ 

linear  combiner  (ie:   tapped-deUy  line)  type  ^  ^ 

more  sophisticated  adaptive  design.  ^™e   c     Qr   recursive 
we  would  use  the  ^T  implanted  with  an  FFT  c  p 
DFT  to  produce  a  series  of  -  P      » -5the  Qut  t  subtracter 
algorithm  would  simply  adapt  the  weigncs  o 
that  would  cancel  the  various  sine  waves. 

»T  'irr^gre^ive  rT.'fhr-ordrof'the 

'^Iter^set  *£"  i  ial  delay  £  -i- ---^ 

equally  spaced  ^^^A  eYu  ut'oflhe  delay  ,-  are  used 

pole  producing  action-  at  the    P       band-pass  filters. 

to  cancel  m  of  these  ,eros  resul   g   ^  narrowband 

The^e  band-pass  filters  in  iu  adaotive  algorithm  uses 

interference  from  the  input  signal.   ™™P^nc(J  in  order  to 

a  gradient  search  to  find  the  correc  «r«  to  ance  1^^^ 

minimize  the  output  power  thus  cancelling 

interference. 
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y(n) 


Figure  1-2. 


Adaptive  Recursive  DFT  FIR  Notch  Fnt 


er 


x  (n) 


Figure  1-3. 


e<n)   (Notch  output) 

y(n)   (Sine  output) 
de(n> 


s(n)  = 


dkj  (n) 


Adaptive  IIR  Signal  Cancel 


er 
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2.0  Superiority  of  the  Second-Order  Bandpass  Canceler 

The  second-order  implementations  of  sec tion  1.  *  «d  *"2  °"er 
considerable  advantage  both  in  ^^h^Sr-order 
adaptive  performance  when  compared  to  ^e  ^£deUy_line  type 

"I'^notch0  ilte  Tu  ch  «  St-  section'  1. 3  would  require 

S^r^-^t  of^on  t.5  would  require  over  100 
weignts  to  accomplish  comparable  performance. 

"Tne  only  higher-order  -tion  that  appearsto  ^competitive 

wlth  a  second-order  ^;      ^Te  required  in  Figure  1- 
approach  is  used,  a  delay  of  n.  made  comparable  to 

2,  but  the  adaptive  part  of  the  *>    t  Df  performance,  the  DFT 

that  of  a  second-order  syste*.   In  term  so  p^    ^       ^ 
approach  of  Figure  1-2  will  create  a  sm 
ripple  in  ^e  resulting  output  sxgny,  but  a^long  ^  ^^ 

laroe,  this  ripple  should  not  be  a  pr  experiments 

that  this  approach  «es.rves  further  attention,  o_    P     ^  ^ 
indicate  that  the  second-order  band  pass  can^*   PP;uture 

best  choice  for  our  P-^  app^UonS;  dements  in  such 
developments  in  rea  -  i    FJNhardwareithmetic  £      ^  ^ 

circuit  of Figure  1-2  might  make  this  approach  attractive. 

2.1   Cascaded  Second-Order  Adaptive  Notch  Filters 

The  filter  of  section  1.1  was  fomented  and  tested  to  see  how 

it  performed  in  the  elimination  o  sin* .waves  in^the^ 

Saussian  noise.   Initial  tests  ™  ine  wave5  are 

approach  performed  very  well   however   w he  . 

close  together,  a  problem  arises.   Fiaur*  -  '        0.125f«,  and 

used.   This  input  consists  of  3  sine  waves  at  0. If,,  B.  ^  .^ 

S.375f,  (36",  45",  and  135")  in  Gaussi an  noise   Ea   sine  «.ve 

and  the  noise  have  the  same  energy  in  the  .ign.l  <^ne  wa 

amplitude  1.414  and  the  noise  is  Gaussian  w  th  variance  l 

Finure  ''-lb  shows  the  output  of  the  first  stage  aT 

f liter  Figure  2-lc  shows  the  output  of  the  second  sag, 

F  gure'2- Id  shows  the  output  of  the  third  stage  of  the  filter. 

Instead  of  eliminating  completely  one  sine  wave,  the  first 
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half  way  between°+h!]  *llter  (shown  in  Figure  2-)M  ,-  . 

of  the  sine-   ?h  V"°  sine  w^es  rather  thin  ,  • Centers  itself 


;-2  The  Kwan  and  Martin  Filter 


Meeting  anS  ZJ^^^^ ,^ ,    the  problem  of 

"s  eama  neateTb  °aeDimn°ate  th*  sinuso  &!*£ £«**   "Iters 

vers:  on  of  one  of  the  Pa55  'lUer  wh°se  output  is  a  '  *inus°i°s 
can  oerform  both  t   ^^H'5-   ^"ce  this^m"^^"^ 
configuration  which  !  5   nafS1"9le  Captive  filter    Structure 
extremely  well.        Sh°m  to  °*  highly  robust  and  perform 

the*  Z^^ZlTs^l   *"2  ""  «-**"  «ructur.  t.  that 
configuration  i.",?^ {»  £»-'"  the  -Pti~  £  eV^ 

sinusoids  to  be  detected  and  r=lo  '„  ^  N   is  the  "umber  of 
>n  rea,-tlme  situations  wm  morTth   T^  beCOmes  ^aotical 
geometric  increase  in  the  reauTrL ^  -  sln^oids  due  to  the 

eou.es  the  required  hardware  to  3N-I  „    -  structure  that 
the  performance  of  the  system.  QUadS  With  no  effect  on 


Ei-i<=> 


EjCz) 


V,  (2) 


^9ure  2-2.   One  Section  of  K 


wan  and  Martin  Notch  Filter. 
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2.3  Kwan  and  Martin  Structure 

The  Kwan  and  Martin  structure  consi sts  o fa  c. scade  of  IIR 
notch  finer,  one  stage  of  which  is  ^^J^Sc  ^  tT 
stace  consists  of  a  bandpass  filter  «th  ^  Su, 

M^-Sr^r^lSirn.l-i^  tranter  function: 


1     2  1  -  -'2 

l~ri         *  (2-1) 


HBP;(Z)  =  7"   -i       2,-2 

BPl        2    1  "  2rico&ei  z    +  ^  ^ 


where: 


r.  =  pole  radius  of  the  i-th  section 

8  ■    -   2ttUj  /u, 

u*j  =  peak  frequency  of  the  i-th  section 

u'  =  sampling  frequency 

-  -  Evsrjr^s-TLS'S'-K  ays 

£  ^rtlntTand^dfh-filter   in  whi, ,  the  pole --££  »  • 

constant  and   only  the  frequency  »,    ^s  f^^^   in  Kwan   and 
approach  which   keeps  a  constant   Q   is  also  d*='"ss  th 

u     ..■  irithpr   of    these  approaches  may  be  used   wicn 

t"ct,,eE:r:nr;opos1n  %uho,?h  r  T  r:  Tsl^ 

concentrate  on  the  constan .   a; ;d« dth  ca «  in  our^  ^ 

addition  Kwan  an  J  "arUn  selec     .^  ^  ^ 

way  -hat  it  is  fairly  ea  >  t   WG  have  not  followed 

the  adaptive  parameter.   For  5imPilC*  J?     narame+er  k-cos^i- 

u   k,,+  haw  rhosen  the  adaptive  pardmewci   j      i 
this  approach,  but  have  cnosen  ™  viPld=  the  following: 

Substituting  this  into  equation  (2-1)  yields  tne 


-2 


1-rj        1  (2-2) 

HRPi(z)  = Ti  2  -2 

EPl        2    1  -  2rtkA  z    +  rx  z 

c    F.nure  2-'  we  see  that  the  notch  filter  for  each  section  is 
the  dH'erence^between  1  and  the  bandpass  filter,  hence: 

(2-3a) 
HHi(z)  =  1  "  hepi(z) 
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Hni<z)    = 


_    1+ri2    1    -  4kiri2-1/(l 


+r:2)  +  z-2 


1  ~  2riki  z-1  + 


ri^ 


-2 


(2-3b) 


CALCULATION  OF  THE  GRADIENT 


The  basic  structure  of  the  k 
shown  in  Fagure  2-3  is  a  Mc_i 
riqure  2-2. 


a  casca 


Kwan  and  Martin  adaptive  fiIt 
de  of  N  sections  of  the  for*  n 


er 

orm  of 


Fi  nnrp  ~>— t    i.' 

y    -  -'■   ' wa"  and  Martin 


Notch  Finer  for  Tracking  N  s:nusoi 


ds 
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The  ever a 


11  transfer  function  is  given  by: 


n  (2-4a) 

T(z)  =."  HNi(z) 

i  =  i 

n  (2-4b) 

=  H  (1  -  HBPi) 
i  =  i 

«h  Martin  choose  as  their  objective  function  J(z)  the 
SSr^ The  oufpuc  o<  the  fln.1  stage  o,  the  cascade: 

J(z>  =CE1(z)]=  =  T(z)2X(z)=  (2"5> 


Hence. 


,  the  gradient  o*  the  objective  function  J(z»  is  given  by: 


&J(z)            dT(z)  (2-6) 
=  2E1(z>X(z) 

respcit  to  each  k-: 

dT<z)   N       aHNj(z)  (2_7) 

=.n  hNi(z)  

From  equation  (2-3a)  we  have: 

anNj(z)   aci  -  HBPj(z)3  ^  -aHBPj(z)  (2_Q) 
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From   equation    (2-2)    we   have: 


D(z) 


-  HBp.  (z)H£    (2) 


(2-9) 


where: 


D(z) 


=    '    ~   *j"j«-«    *   r. 


He:    (2)      = 


Cd 


©nominator    o-f    H 


BPj  <2>3 


Substituting   equation    (2-9)    int 


°  ^uati°n    (2-7)    we  ODtain: 


dT(z)      N 


•  Fj<2>HSj(z) 


(2-10) 


notch  Alters  without  the  notch  f Iff    Sener.ting  the  product  of 
equation  (2-10).   To  gene"t=  th         J  '  Whlch  is  required  in 
would  require  N-,  biq'abTpIr  se  tion^V^  ^h  ^Uon 
M'N-1)  fai quads  just  to  generate  fh    resu'ting  in  a  total  of 
able  to  reduce  this  by  u^ing  th.   VT^*"   ^  and  M*rti"  -re 
the  input  to  their  cascade?'  Since  th^  V^  bandpaS5  «»- 
of  the  required  H  .  (-)  *=J        thls  outP<-it  already  has  (i-ti 
Squads  are  needeS'for  1^1   £"'  °^    W^    additional 

N 


^  (N  -  j)  = 


N  - 


j=i 


j  =  i 


=  N2  - 


N2+N   N2-N 


(2-11) 
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Addir.q 


this  to  the  N  biquads  required  to  realize  the  cascade  of 


Iters  and  the  N  biquads  require 
"actors  yields  a  total  number  of  biquads  given  by: 


d  to  realize  the  H5j (z) 


N2-N 


N(N+3) 


#  o-f  biquads  (Kwan/Martin)  - 


+  2N  = 


(2-12) 


2.4  Improved  Structure 

Fipure  2-4  shows  the  improved  adaptive  notch  ^^er  structure 
l-igure  -  ">  „   Tho  i.ev  to  the  improvement  is  the 

proposed  in  this  paper   T^_k^.I°(,^(z)  for  the  cascade  of  the 
recognition  that  the  output  E,  <;>-^X<;;   t  of  the  individual 
notch  filters  can  be  written  b°th*s*  P™,  time  the  input  and 
notch  finer  section  transfer  functions  HN<->  time       ^ 
in  terms  of  the  input  X(z)  minus  the  outputs  YjU)  of 
bandpass  -filters: 


H51(z) 


HS2(z) 


4 


T(l 


-H 


d; 


-H 


H^z) 


-H 


H5N(z) 


-> 


Figure  2-4, 


New  Notch  Filter  for  Tracking  N  sinusoids 


-ie- 


Ei(Z)  =T(2)X->  -J.H^C.)  X(z) 

1-1  <2-13a) 


X  ( z  )  -   ')  y .  / ,  i 

ilt  <2-13b) 


To  get  the  product  of  H   (?i   ■* 


we  may  use 


i!iHNi(2j  X(z)  =  T(z)XCz)  -  Y  (z) 

J  "'  (2-14) 


**4 


Figure  2-4  makes  use  of  th ■ 

e"tr*  HEPi(2,.   Thus  we  have:        '^t  stage  does  not  need  the 

*  °f  biquad£  "*  structure)  =  3N-, 

(2-15) 
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Table  2-1 


2.5  Experimental  Results  Using  the  New  Structure 

The  algorithm  was  compared  -  the  K«an  and  .arUn^lgor.thm 
and  found  to  perform  as  wel  or  better^ *  ^   ^  ^ 

all  of  the  experiments.   Figure  ~  -  waveB 

algonthm  with  5  notches  applied  to  a  signal  ,xt  ^^ 
imbedded  in  noise.  Figure  2  oa  snows  Fmterinq  and  Figur<= 
response,  Figure  2-5b  shows  the  output  «f  .««.  filtering. 
2-5c  shows  the  frequency  response  °<t*%°fj3"lrl2ters.  Due  to 
Figure  2-5d  shows  the  adaptation  of  the  five  P™£ 
computational  complexity,  *  ^^  D°el  °  ^  to  implement 
the  Kwan  and  Martin  algorithm,  thus  just  oe  g  Martin. 

the  new  algorithm  represents  an  improvement  over  Kwan 
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«"£^2  JET^t  i«r-  -«"  «-  5-notch  adaptive 
Figure  2-6a  shows  the  input  spectrum  p""  lmbedded  in  noise, 
output,  Figure  2.6c     P  th^^r-J*""''-  2-6b  shows  the 
shows  the  adaptat1Dn  of  the  5  notches^  r,   ''  ^  Flgure  2"^ 
cancel  the  three  sines,  while  the  rll' ■  **  °<   the  notches 

around.  '  wnile  the  remaining  notches  wander 

mteFr9isreus2ed  wheTtnll  ST!.*" .«»  *"»tch  adaptive 

2-7a  shows  the  input  spectrum,  Fi aTe  ?%"?   '"  —   F>9-e 
Figure  2-7c  shows  the  output  «n=r?    ""   Shows  the  output, 
adaptation  of  the  5  notches   Shat  happen"  Fi°"re  2"7d  =now5  the 
interesting  and  shows  just  how  robul^K    J"  FlQUre  2"7d  is  very 

^irst  4  notches  lock  solidlv  on  th=,f  fl90rithm  is"   T"e    * 
K°,    60°,  and  90°).   The   L  n  ?'"  "  '^'endes  (30°, 

refining  three  frequencies   2S°   ns°  '!"  ^ln»  "-tn.en'the 
unable  to  cancel  all  three,  1 1  spends  .'if??,  T' '   Since  4t  ls 
each  one,  thus  attenuating  each  o th=    n   *  5lt  0f  time  °" 
scan  across  the  entire  range  in  hope  of '  .0ccasi°nally,  it  will 
solution.   But  it  quickly  lettlJl T  ?f/  ndln9  a  better 
attenuating  each  of  the  Gaining  sfnes'°  ^  Patter"  °< 


"g^roV^^He^r^^eo^^f"1^  the  Slne  ««-  i-  snow 
demonstrate  the  abi  fityl^^^  "  sines  in  order  to 
waves.   Figure  2-8  shows  the  Ln?  t     ^°  process  10  sine 

the  adaptation  of  the  10  parameters. 


2.6  Conclusion 


We  have  demonstrated  that  fkn 
structure  uses  considerable  ess  h^H    '"^  n°tch  <ilter 
proposed  by  Kwan  and  Mart  n  vet  nf,   "  *"*"  the  structure 
typical  applications.   In  addition  ,1^  ^  Mme  <0r  bet^r) 
here,  chapters  3  and  4  of  thl5  1   *"!  reSuIts  demonstrated 
Provide  other  experiments  to  verffv  t\    *"  ""^  TheS1S  CW 
structure.   This  work  has  L  w  V     Performance  of  the  new 
"0.113.   In  all  of  these  cases  fw  reported  at  two  conferences 
-11  or  better  than  tnTL^d  ^st^ref  >*""-   " 


in 
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Input    (5  sines   in  noise) 

— t 1 — ■ 1— 1     ~" 


.— -_.  -riy.v. 


"~ 1 


I^A-»^/^./v^AV 


.2         .25         .3 
FFEQUENCY 


.35 


.4 


.45        .5 


J 


Figure   2-5a      Input    o*    5   S.nes    (30°,    45°,    36°,    60%    and    135°)    in   Noise 
Number   of    sinusoidal    inputs   is     S        A«F    -    a. 414,    1.414        ^        ^  ^    ^    . 
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-^^ii^lELAdaptiue  Filter 


U 
a 
I 
u 

e 


Figure   2-Sb      Fi  ]  f  n 

Input   3  Sines    &?  *£put°f   5-Notch  Adaptive  Filter 

»-■•»«..      slimit      B  ' 

tmg.      Factor   =    {0.900,    0.9B0i    ^^    ^    ^ 
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Output  5-Notch  Adaptive  Filter   (5  sines    in  noise) 


^V-'.^'V 


.../-••  Vv/V 


V  V 
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.05 


15 


,2         .25  .3 

FREQUENCY 


,35 


.45 


«-  ^T^^r^.^Ksw^8^ FUte 

Number  of  sinusoidal  input.  «  -»   «*  "  '1'4^  i"  {0.J  0.0,  0.0,  0.0 
FREQ  =  C.B83,  .125,  .167,  .*j0,  .*'<*'  =   1,    n  sec  i 

Second  order  forgetting.   Factor  -  ,0.900,  B.viflio,  ».™  i 
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-^Ill^L^LI^^eters  for  5. 


Notch 


a 
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11 
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Iterations 


_OC_ 


Input   (3  sines    in  noise) 
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358    - 


0  .65 


.1        .15 


.2         .25  .3 

FREQUENCY 


,   „.  ,-rao      -?a°      ancj    135°)    in   Noise 

Figure  2-6a      Input   o-f    3  Sines    (.0    ,    ^6    , 

limit    -  0  000000*        <™in   =0-000001?        pmaX    =    1000e 
step   =   0.100000*         5-hml      -.0"BT   ILp   -a    414,    1.414,    1-414} 
Number   o-f    sinusoidal    ^tslt  ^      ANG  .    {0;000,    0.000,    0.00B> 

FREQ  =    (0.083,    0.125,    B-o75>        ANb  =  n_sec   = 

DC   .   0   000  NOISE   -    I    «■    ,  - ^^1.    00,000i    0>000} 

Second   order    -forgetting.      Factor   -    J8.WU,    *■         < 
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Output  of  5-Wn^flJ.r,ivrr||[||    (3 


nes    in  noise) 
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we     158     2ee 


He^cf"       35B       «       «       580 


Figure   2-6b      Pi  l  f 
Input   3  Slne5      ,00        °otPL,t   °f   5~^tch   Filter 

-fo  n  "  •'   and    135°>    m   Noise 

»tep   =  0.130000  . 

otch  pole  -dlr=k  0;9^-00009000-t^-b0b'  ■■«».  « "«■>    -sec  =  5 

*ond  order  forgetting   V.rf'  ,1      '  0-'00'  0-900J 

-    ng.   Factor  =  {0.900i  0_900j  0_900i  0900]  0  90^ 
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Output  of  5-Notcn  Adaptive  Filter   (3  sines    in  noise) 
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FREQUENCV 
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Fiqure   2-6c      Frequency   Response   oi    Filtered   Output   o*    5-Notch   Filter 
Input    3   Sines    (30°,    36°,    and    135°)    in   Noise 

step   =   0.100000,        sjimit   =  0.000000,        pmin   =0.000001,        pmax    =    100000. 
Number    o-f    sinusoidal    inputs   is      3:      AMP   =    a. 414,    1.414      1.414, 

FREQ   =    C0.083,    0.125,    0.375}         ANG   =    C0.000,    0.000,    B-BW>  _ 

DC  =  0.000,        NOISE   =    1.000,         kjimit   =    1 .  000,        slow  =      1 ,        n.sec 

^dVtg^rVS  e-900,    B-900.    »-™.    ■■"»> 
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-^m^-r~ZZ~^        r-J^__^il£L_C3_SMiesj in  noise) 


V 

a 

1 

(i 

e 


8     "    '"    '^Z™J^-^-^-15    m 

figure  2-6d      p3^ 

00      parameter   AHa„4.    ±.  ■ 
input   3  Slnes    (30o,t?   ^-  ^    S-Notch   Fnter 

eP   =  0.103000  c    ,.     .  13"'    >    ln   N°ise 

nber   of   s-ni     ',,    f-fimit   =  0.000000  . 

ppcn   S-nusoidal    inputs    k      ?        ?'        pmin   =0.000001 

[*i*i  vaues    s,V-,r„'    kji.?r- 1 10000,  0;000'  0-000 

Kpoie  ,adir=ki0;90  ;T^-T^-0  v:000  ^W-    n-5ec  =  5 

-ond   order    forgettino       Vlr*  "'    0"900'    0-9B0> 
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Input  7  sines    ( 


(30 


,   45,   68,  98,   128,    135,   &  156  deg)    in  noise 


708   .. 


680    ■- 
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K^^Aa-^^^  wJL^  ww* 
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.2         .25       .3 
FREQUENCY 

Figure   2-7a      Input   o-f    7  Sines 
(30%    45°,    60°,    90°,    120°,    135°,    and    150°)    in   Noise 

step   =  0.100000,        s.limit   =  0.000000,  pmm   =0.000001,        pmax   =   100001* 
Number    o-f    sinusoidal    inputs   is      7. 

AMP   =    (1.414,    1.414,    1.414,    1.414,  1.414,    1.414,    1.4   4, 

FRED   =    (0.083,    0.125,    0.167,    0.250,  0.33:.,    0.375,    0.417/ 

ANG  =    (0.000,    0.000,    0-000<    0-000'    0*000'    0'000'    0'00 
DC  =  0.000,        NOISE  =    1.000,      k.limit   =    1.000,        slow  =      1,        n_sec 
Initial    values   o-f    kl    =    (0.000,    0.000,    0.000,    0-000<    0;000J- 
IMntch   nol^   radii    =    (0.900,    0.900,    0.900,    0.900,    0.900) 
Secor/orcerlorgetting.    Vactor'=    (0.900,    0.900,    0.900,    0.900,    0.900 
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^^XI^£}L^Pti^  Filter   (7  . 


LULJZWt    in  noise) 
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'\™Uo™      ^~^~^~^ 


p  ■  •  "-c  «  l  iuns 

Input  7  "^^Vi^^g.^*  «»t 

5  =  0-100C00,        s   ,,    ..        n  *"      '    i-J   '    and   150°)    in   Noi 

^  °f,rnU501d^    -   uts'irr-        Pmi"   =°«1  on 

AMF    =    n.414,    1.414,    1.4i4S    /l,  .  '        DmaK    =    1000000.00 

FREE!   =    {0.08,-t,    0    ,,t'    „  ?;,'    1-414,    1.414.    1    41.      ,    .,„, 
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Output  5-notch  Adaptive  Filter   (7  sine    input    in  poise) 
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FREQUENCY 

Figure   2-7c      Frequency   Response   oi    Filtered   Output   oi    5-Notch  FiltJ 

lit   7   Sines    (30%    45%    60%    90%    120%    135%    and    150°)    in   No.se  . 

step   =  0.100000,        B.limit   =  0.000000,  P-in   =0.000001,        P«™    =    100000 
Number    o-f    sinusoidal    inputs    is      7. 

AMP  =  CI. 414,  1.414,  1.414,  1.414,  1.414,  1.414,  1.4  4) 

FRED  =  '0.083,  0.125,  0.167,  0.250,  0.333,  0.375,  0.41./ 

ANG  =  C0.000,  0.000,  0-000,  B-000<  0'000'  0-000'  0'T} 
DC  =  0.000,   NOISE  =  1.000,   kjimit  =  1.000,   slow  -  1, 

Initial  values  ot  kl  =  <0.000,  0-000,  «-000'  0'000'  0;000> 

Notch  pole  radii  =  C0.900,  0.900,  0.900,  0.900,  0.900/ 

Second  order  forgetting.   Factor  =  -C0.900,  0.900,  0.900,  0.900,  0.900, 


n  sec  = 


-TO 

'Ji. 


IJT>,      r,       r**n    ,  r-H-V811*  inputiTi  noise) 
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106       268 


460         568       668 
Iterations 


788      see      gee     lm 


'         ^Put   /&£  g?§   £-£r  Adaptation  ^5-Notch  rate, 

^P   -   0.100000  _    ,,     .,  '      '      '      ^    '    and    I50O>    *n   No1Se 

Mumber   of   sinusoidal  \nn   T  0"000000,        pmin   =0.000001 

AMP  =    ^    414      7   i,inputs   1£      7.  .000001,        pma.,    =    1000000_00 

freq  =  ;0w  0-^4>  i-«m.  i.4i4,  i.4i4,  1.414  ,  4M, 

■  ~-:*Kiv3. r  •  ••?•'  '•-••  Si ; £ 


Overlay  ofPara*et^^ 
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288         258       380 
Iterations 

Fiaure  2-8      Input   o-f    10  Sines  in  Noise 
<30%    36%    45%    60%    72°,    108%    120%    135%    144%    150°) 

5teo=0. 100000.        s   H^t    =   0.000000,        P-n   =0.000001, 

"£"-  °1    JSwS/l.W,    1-414,    1-414,    1-414,    1.414      1.414,    J 

DC  =   0.000,        NOISE   =   0.000,        k_li»lt   -    1.000,        slow 

ln%*lw.TmZ  0!000,  ■■«»,  °-000-  0-000'  0-000'  0-000'  B-00B'  0'0B1 

NDtf0  P900:   O?d.00,=0.900,    0.900,    0-900,    0-900,    0-900,    B.™,    0-900'    0"W 

^V^^IS  ££%%«.  ■.*». ..«,  ■.'«. 0-^ 0-90> 
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3.0  Detailed  Description  of  the  New  «goritnm 

In    this   chanter    n-F    fi- 
new  algorlthm  and  delstratT"     "J"!  d™£   in   "*•»    the 
The   derivation   of    the   structure   and   fn*   °D   simul^ed   data. 
discussed  in  sections  2.3  and  *  4       Mp"  9radient   "Elation  were 
details   for   fomenting   the  ad^at^nn"  "V^"    9We  the 
choose  the  various  adaptive  otrTj  '    lnf°™ati°n   on  how  to 

effects  of   different   cho  ces  on   the  oe^f"'  dem0nstra«°n  of   the 

'   tne  Performance  of   the  filters. 


3.1    The  Adaptation  Algorith 


m 


^ucFtu^co^ta?nsVetS„ohtypSesrUoftU:!  <°:  ^  -  al9orithm, 

types   o-f    second-order   filters: 


The 


Band-pass   filters    (see 


'quation  2-2)  : 


H..:  (z)  = 


(3-1) 


and  Sensitivity  filters  ( 


see  equation  2-9) 


HSj(c>    = 


£T     2 

J 


-i 


1    -   2r.k.z"1    + 


J    j 


r  . 


(3-2) 


where : 


b(z) 


=    1-2 


J    J  j  '■deno.7!i,p}af  "»-   0*   u 


'?>  1 


The  partial  derivative  is  given  by  (see  equation  2-10): 


dT  (Z  )         kj  /T_7\ 

=  -.1  HNi(z>  HBPj(z)H5j(z)  ,3S 

-  1  i  =  1 


The  only  parameter  that  is  modified  in  the  adaptation  procesl 
is  the  pole' cosine  k-  -for  each  bandpass  -filter  section  J.   This 
parameter  also  effects  the  denominator  of  the  sensitivity  filter 
H  .(z>.   The  only  other  parameter,  r-,  is  fixed.   Selection  of  rj 
Hill  be  covered  in  section  3.3.   The'standard  LMS  algorithm  would 
give  us  the  following  update  procedure  for  adapting  k  ■: 

dT(z) 
k:.=k  -M.(z) (3'4) 

V 


w 


here:  u   is  a  fixed  step  size  (typically  #=0.1) 


As  derived  in  chapter  2,  the  correct  choice  for  the  error 
function  e(z)  is: 

=  (Z)  =  T(z)     (for  all  j)  (3"5> 

j 

However,  an  alternative  choice  which  can  be  derived  by  forming 
the  gradient  of  the  section  outputs  YMz)  is: 

ej(2)=Vj(Z» 

Both  of  these  choices  for  e^z)  were  tested  and  found  to  work 
well . 

A  word  of  explanation  is  in  order  here.   Equation  (3-4)  seem 
to  imply  that  we  are  mixing  z-domain  and  time  domain  terms.   in 
sens.,  we  are.   The  derivative  are  taken  on  the  frequency  domain 
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SSl^ti^^^'  |f  ?«  "<*>'">    in  Figure  2-4  are 
equation  <3-4>.   A]sD)  S.1^"-^-.  series  that  a,^  in 

Vj(i).   This,  of  course,  imnlies  thi/     really  the  time  series 
values  to  represent  the'a      a  f  "f  "re  """9  instantaneous 
approbations  to  the  actual  values!  '    *  these  are 


3.2   Forgetting  Fact 


or 


*nff  SSt    Vari0L,S    ^rqettinq    L^"         ^    (<?g:    See    C5J).       In 

Allowing   update   equation:  **=torS|    WG   employed   the 

^n::]:=kjEk,_5tep 

ll-w+P^»03/pmax)/(pmin+prk?0;i) 


where: 


step 


=  *  (fixed  for  a  particular  structure) 


(3- 7a) 
C3-7b) 


pCk,0] 


=  output  of  forgetting  filter 

(•-■-7c) 

=  limits  the  effect  n*  7 = 

ettect  of  large  pCk,0]        (3_7d) 


pmin 


-  limits  the  effprf  ^x 

G++ect  of  small  prk,03 


<3-7e) 


NOTE: 


We  found  that  nma»  anpw 
«-  ^Pnate  algori t^ouTc!  FJ^V**™-  "  '"  "«"* 
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k j C  k ] : =k  j [ k ] -step*gr ad i  ent/pE  k , 0] 


(3-8) 


However,  all  of  the  data  in  this  chapter  was  ran  with  equation  7 
setting  pmax=1000000  and  pmin=0. 000001 . 

Three  different  forgetting  filters  were  tested.   The  first, 
called  zero-order  forgetting,  uses  the  sum  of  the  squares  of  the 

stantaneous  outputs  »TCz)/akJ  as  pEk.BJ.   The  next  -Ued 
first  order  forgetting,  is  a  simple  «n9le-pole  first  order 
DaSc  filter  applied  to  the  eection  output  CATD/Ak^  .   m« 
third  type,  called  second-order  forgetting,  is  a  second-order 


low-pass  filter  applied  to  CdT(z)/dkj3 


Zero  Order  Forgetting 


pCk«0]  = 


(3-9a) 


First  Order  Forgetting 


pLk,0J  =  d-rf) 


dT(z) 


ak 


+  rfpLk,03 


(3-9b) 


Second  Order  Forgetting 


pLk.,0]  = 


output  of  one  of  the  following  second-order  filters 
with  input  CdT  (z)  /dk  ■  32: 


-38- 


Simple  2nd-0rder   Low-P 


ass 


<3-9c) 


Complex  2nd-0rder  Low-Pass 


h,(l  - 


H2(z)  = 


where:       S.    =   ArcCosO, 


OS<5£ . )2_1     +     ,-2 


(3-9d) 


The   filter    oi    equation    CS-Orfi 

3-3  Selection   D?  **£  Bandwidth5   ""^  ^  ^   and   ^«n    T53. 

Each   bandpass   filter   h 
parameter   r      '     H   we   let   ,  bVtnTf*  ban?width   ^   by  the 
notch   and   yJbe  the  upper   Mb   «™        "^   ^   fre^ncy  of   the 
bandwidth    i«   d  =  ,,    _ 'P;     wfshf:f?U?nCy  of   t^  notch,    then   the 
frequency   z=j    <9sc         *      We  ^11    flrst    calculate 

<*her    vaU,es  o,    ./ba^'on   £", ^f^    «*«*.*  .   and  y  for 

^"S  ^x^^r the  tw° 3db  *-q«es  K.  and  „ 


=  |HN(z)  p  = 


1+r 


2  12l 


d+cos2ut)2   +   sin2 


2ut 


<l+rWUt>2   +  rW^7'        <3"10a> 


1       = 


4   +   4cos2o)t 
(1    +   2r2cos2yt    +   r4 


(3-10b) 
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^.  2^  * 


(3-1D 


cross-multiplying  yields: 

1+2r*cos2.t+r«  =  l+2r2+r<  +  c«artd*V) 

U.ing   the   fact   that   cos2.t   -   W-t-1,    we   solve   equation    (3-11) 

■for    cosut: 


cos2ut   =   2cos 


2ut   -    1 


-1?* 


1+r 


(3-12a) 


1-r 


(3-12b) 


cosyt 


=  + 


-\  2U+r4) 


Since  xB   is  the   lower   3db   .regency  and   y.   is  the  upper   3db 
■frequency,    we  have: 


=   ArcCos 


-(1-r2)     1 


JiTl+r* 


(3-13a) 


r       1  — 2     "» 


y      =   ArcCos 


-J2(l+r4 


(3-13b) 


Also 


■    2fl    .    rn<z29   =    1 ,    we  have 
,ince   sin   o    +   cos   u         *  ?    " 


1-r 


cos  (>■:«)    =   ~ 


2(l+r4) 


(3-14a) 
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COS(y-)  =  + 


1-r 


2(l+r4) 


(3~14b) 


sin(x0) 


=  + 


1+r 


2(l+r4) 


(3-14c) 


sin(y0)  =  + 


1+r 


2(l+r4) 


(3-14d) 


WS  "*  h°W  Mlcu1*-  «»  bandWIdth  d  from  d  = 

COS(d)  =  CDS(Ke-ye) 

=  CO.(x.)co.(y.,  +  .i„,x,,Bin(yt, 

Substitution  fi 


xe  ~  V.: 


rom   equations  (3-14)  int 


(3-15) 


°  equatlon  (3-15)  yields: 


cos(d)  = 


-n-2r2+r*)  +  (i+2r2+ 


r4) 


2r2 


2(l+r*) 


1+r 


(3-16) 
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[l+r2l 


r  1+r 


2  1 


-  2r  k  •  2  1  + 


L  2     J 


z"2  =  0 


(3-17a) 


i+r 


L  *    J 


4rk  ■ 


1  - 


1+r 


i  .-1  +  *"2 


=  e 


fT-1Th) 


Substituting  ^or  z=e<i f =cos </r ) +  jsin (/r )  in  equation  <3-17b)  and 

setting  the  real  and  the  imaginary  part  to  zero  yields: 


1  - 


4rk  ; 


1+r 


cost/r)  +  cos(2/r)  =  0 


4rk  . 


1+r 


sin(/r)  +  sin(2/r)  =  0 


Using  trigonometric  identities: 


(3-lBa) 


(3-18b) 


4rk  : 

1 1  cos(/'7)  +  2cos2(/r)  -1  =  0 

1+r2 


<3-19a) 


4rk 


—  s 


in(/r)  +  2sin(/r)cos(/r)  =  0 


1+r' 


both  equations  yield  the  same  result: 


:os(/7)  = 


2rk 


1+r 


(3-19b) 


(3-20) 
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^•.^^"^^^^^"ly  -lightly  different  ,rom 

frequencies  would  be  identical  at   r'  r=1'   the  twD 

be  less  than  one.   Table  3-1  iimt°   tK   Sf '  f°r  stahility,  r  must 
vaue  of  r.   From  the  table  we  See  that'!  fti0nShiP  °<  »   to  the 
well  suited  to  the  type  of  notrhff    ,      alues  of  r  ar°""d  .9  are 
b-,d  interference.   fhPe  t°Ibl Tl £  ?^f  *°  eI^-te  narrow-"8 
the  actual  notch  frequency  /  t   ,  „U,  he  measured  values  for 
upper  3db  frequency  J  for  "each  of   he  not  h  '^^  "«  *nd  th* 
closely  to  the  calculated  values  in  „   >  ^   These  c°mP^e 
20',   and  as  follows:  '"  ^t.on  (3-16),  equation  (3- 


Table  3-1 


Effect  of 


r  on  Notch  Width 
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s,/   —  x  =  d 


(3-21a) 


rT-91h) 


hence: 


(3-22a) 


y2  _  xy  _dy  =  y2  -  -f2  -dy  =  0 

y  «  d/2  +  4^^T7     x  =  d/2  -  4(d/2)»W2  C3-22b) 


;.4  Minimum  and  Maximum  Frequency  of  Notch 


Equation  (3-20)  implies  that  the  notch  frequency  is  slightly 
different  from  the  pole  frequency  6.       In  particular: 


)s(/r)  = 


2rkj 


2r 


cos(S) 


(3-20) 


1+r2    1+r 


This  difference  becomes  smaller  as  r  approaches  one   However ,  if 
the  Pole  frequency  5  is  at  its  minimum  o+  0  ,  then  kj-i  and. 


\  =   f 


77 


2r 


=  Arc Cos 


#  =  e 


1+r2 


(3-21a) 


If  6    is  at  its  maximum  of  tt  (180°),  then  k^-1  and: 


tt  -  A  =  fr 


-   ArcCos 

0=TT 


-2r 


.1+r2. 


<3-21b> 


Fro.  .«,.«.„  (3-2BI  ..  ««  th.1  U»  '"'""".TT  V  ""' 
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frequencies  goes  to  zero, 
various  values  o-f  r. 


Table  3-1 I 


gives  A  and  180°-£  f 


or 


Table  3-11   Effect  of  r 


on   A 


r 

"iiiimum  Notch 

"frequency 

A 

Maximum  Notch 
frequency 
180°-A 

0.3   1 

56.60151° 

123.39849° 

0.5   | 

36.86990° 
20.01596°    , 

143.13010° 

0.  7 

159.98404° 

0.8 

12.68038°    1 

167.31962° 

0.9 

6.02558°    J 

173.97442°   1 

0.95 

2.93760° 

177.06240° 

than  the  minima*  notch   frequency  or   ST       *   *  fre^^y  lower 
notch   will    be  -forced   to  this  2imZ     1        ^   1E  e  DC  val"e,    the 
there   is  a  sinusoid  above  th  frequency.      Similarly     if 

frequency,    the  notchtuf  IT itllTVl^^  ^   *   th*  Cist 
Figu,-E   3-1    shows  the  result   of   a  DC  ■?        J"  ma*lmum  frequency. 
process.      The  notch   is  forced  to  6  2^™' .  °"   ««  Captation 
the  radius  of    the  notch   filter   is   rh.  '"   thls  case  b^ause 

Figure  3-2  shows   a  siqnal    at  th«  M  ?   aS  r=0"9-      Similarly, 

notch   to  the  ma*imum  fluency  lA™^^  <°^  & 


3.5  Aig0ritnfn  Perf 


ormance 


the 


Critical  to  the  perfor 
various  parameters: 


/r.ance  of  the  al  qorith 


im 


is  the  choice  of 
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1.  Number  of  notch  -filters, 

2.  Notch  width  r  ■  , 

3.  Zero,  -First,  or  second-order  forgetting  (or  none), 

4.  Step  size  /i, 

3.   Limit  on  the  parameter  k,  and 
6.   Limit  on  the  parameter  p. 

thp  first  three  choices  are  somewhat  independent  of  each  other. 
But  t  he  cho ice  of  the  step  size  *   is  highly  dependent  on  the 
I-    Lop  rhnice^   We  found  that  limiting  k  such  that   l^k.-i 
first  three  choices.   w*  T        However,  we  found  that  limits 
was  very  important  in  all  cases.   nowtver  , 

on  p  were  not  very  useful. 

The  choice  of  the  number  of  notch  filters  has  a  profound 
eff, -   o  Ihe  algorithm  performance,  but  is  "^^J^  * 
our  control.   We  assume  that  the  number  of  notches  is  fixed 
some  number  by.  the  application.   Thus  in  evaluating  the 
algorithm,  we  look  at  three  cases: 


Case  1:  Performance  with  only  one  notch, 
2:  Performance  with  five  notches, 


Case 


Case 


Performance  with  ten  notches. 


alqori thm. 


Table  3_nI  summarizes  experiments  with  the  new 
For  Ihese  experiments,  the  following  10  frequencies  were  used. 


4 


8 

9 

10 


f ■   =  0.08333  f 


i  n 
i  n 


=  0.1  f 


f .   =  0.125  f 


i  n 


f  .   =  0.1667  f 


i  n 


=  0.2  f5 
=  0.25  f 

;„  =  0-3  u 


i  n 


l  n 


l  n 


i  n 


i  n 


l  n 


=  0.3333  f 
=  0.375  f5 
=  0.4  f« 


(30°) 

(36°) 

(45°) 

(60°) 

(72°) 

(90°) 

(103°) 

(120°) 

(135°) 

(144°) 


^n  cases  where  only  one  frequency  was  selected,  the  results 
are  based  on  the  results  with  each  ^^^^ 
8,  12)  taken  separately.   In  cases  wnere  tnr  «r    h 
used,  the  following  five  combinations  were  used  [(1,2,61, 
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encies  are 
6 


-:  •  _< 


Conclusions 


Table  3-IV  shows  the  results  of  +h 
these  results,  we  conclude  that  a  L   ff  ^P^^ents.   From 
multi-notch  filters  and  t  a   '  r  ?ett,nS  ***«»-  is  needed  fQr 
does  not  see,  to  inprove  the  perform    ^  f getting  factor 
factor  which  seems  to  wor(,  qu^te  sat? T  T"  *  '^-order 
order  forgetting-f actor  is  the  nnf!       °riUy-   The  «r«t- 
"Jump"  convergence.   Here  we  refe  V""  ^  ^  not  B*hl"t 
situation  when  there  are  more  i„tlr«    JUmP  conv^gence  as  the 
last  notch  in  the  cascade  °7u^  Zt^5  ?"*   n°tChes  *»«   the 
frequencies  to  be  eliminated   T*h,   *"eB"  the  r«"*lnina 
actually  be  desireable  "n  some  £*'>  "."mergence;  „,iqht 
«*  to  avoid  the  firSt-orderforngSXS'fLntoWriCh  "**  ™  ""» 

P»ax  m  order  to  obtain  I         specific  values  for  n      °! 
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• 

• 

Table 

3-IV 

Summary  o-f  Tests 

on  New  Algorithm 

■ 

Case  ( 

#  o-f 
2-nd 
Drder  ! 
sezs 

Limit  on 
p  value 

1 

No     Zero-Order   lst-Order  2nd-0rd- 
-orgetting  Forgetting  Forgetting  Forgettif 

P.in  1 

3«ax 

U 

iter 

a 

Lter 

il        : 

Lter 

il       itj 

la 

1 

1 

0.05 

20 

NO 

.02 

50 

.05 

67 

.05 

20 

.05 

Ian 

1 

1 

0.05 

20 

YES 

.02 

50 

.05 

Jump 

.05 

20 

.05   14 

lb 

1 

3 

0.05 

20 

NO 

.02 

700 

.05 

Jump 

.05 

30 

.05 

lbn 

1 

3 

0.05 

20 

YES 

.02 

750 

.05 

Jump 

.05 

50 

.  05   :K 

2a 

5 

3 

0.  10 

10 

NO 

1.46 

2500 

.017 

340 

.063 

200 

.07   <N 

2an 

5 

3 

0.  10 

10 

YES 

1.46 

150 

.017 

1200 

.063 

400 

.07 

2b 

5 

5 

0.10 

10 

NO 

1.46 

325 

.017 

1410 

.063 

150 

.07   1 

2b  n 

5 

5 

0.10 

10 

YES 

1.46 

3120 

.017 

3500 

.063 

175 

.07 

2c 

5 

7 

0.10 

10 

NO 

1.46 

Jump 

.017 

1667 

.063 

450 

.07  Jffl 

2cn 

5 

7 

0.10 

10 

YES 

1.46 

Jump 

.017 

700 

.063 

167 

.07  Jl 

3a 

10 

7 

0.02 

50 

NO 

1.46 

2000 

.02 

1667 

.075 

500 

.07  II 

3an 

10 

7 

0.02 

50 

YES 

1.46 

9500 

.02 

5000 

.075 

750 

.07  | 

3b 

10 

10 

0.02 

50 

NO 

1.46 

7500 

.02 

7500 

1  .075 

2000 

.  07  !5 

3bn 

10 

10 

0.02 

50 

YES 

1.46 

100C 

.02 

1002 

.075 

200C 

1  .  07  :>i 

1.  In  all  ca.es  without  na.se,  parameters  converge  to  exact  value  when  n 
of  notches  is  creater  or  equal  to  number  of  sines.   Iterations  listed, 
table  represent  convergence  to  three  decimal  places  of  correct  value. 

2.  In  all  cases  with  no.se,  parameters  oscil late  around  exact  v,luj.   J« 
number  of  iterations  listed  in  the  table  is  the  number  of  iterations 
required  to  establish  this  oscillatory  pattern. 
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"  ^^r^Z^r^^r   £"  tha"  "^  "-  notches 
table  by  the  word  ..Jump,._        wav«.   Thls  pattern  is  indicated  in  the 

4-   Values  o-f  D  •   anH  « 

«hlch  generally  ",f'"'n  the  tab^  refer  only  to  zero-or*   , 

P.in-10   and  P-ttx=1046)B 
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4.0  Algorithm  Performance  and  Simulation 

The  adaptive  digital  filter  algorithm  as  Ascribed  in  the 
earUer  chapters  is  simulated  and  tested  using  synthetic  data. 
s!m,  ation  was  carried  out  on  VAX  11/785  computer  using  Fortran 
77   The  synthetic  data  generated  for  testing  this  algorithm  is 

o-f  four  categories: 

i.  Sinusoidal  siqnals  with  white  gaussian  noise 

?.  Narrow  Band  Noise  with  white  gaussian  noise 

V  Bi-Phase  Shi-ft  Keying  (BPSK)  sequence 

4*  Frequency  Shift  Keying  (FSK)  sequence 


4.1  Sinusoidal  siqnals 

For  qeneratinq  sinusoidal  signals  placed  at  different 
normalized  frequencies  e,  a  second  order  AR  process  given  as. 

wa=  used  with  Poles  on  the  unit  circle.   This  approach  was  chosen 
in  order  to  reduce  the  computational  burden.   Using  equation  4-1, 
we  find  that  initial  conditions  are  very  important  and  they  are 
chosen  such  that  x  t=0  and  x_„=-sin (0 t >  giving  a  unit  amplitude 
sinusoidal  siqnal/Vhe  B,    value  is  between  0  to  180  and  n  is  the 
number  of  frequencies  desired.   The  required  signal  yk  needed  to 
input  into  the  adaptive  algorithm  is  given  as: 


i  =n 

-  V  v1  +  r  (4_2) 

i  =  i 


where  ru  is  a  white  gaussian  noise  N(0,<r  ) 
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4.2  Narrow  Band  Noi 


se 


The  narrow  band  sinna]  ic  «« 
equation:  lgnal  ls  generated  using  the  difference 

2     *  (4-3) 


The 


i  =  r. 


i  =  i 


(4-4) 


4.3 


Bi-Phase  Shift  Keying  Sequence 
The  generation  of  BPSK  cinn3i 

are  generated  as  shown  "°£ ur^-l^rf  V^*  P-ts  which 
follows:  yure  H    1-      The  three  D^fc  ^„  . 


ee  parts  are  as 


li-iter.   (An  importanyt  notVift^he  -T^  *   hard 
the  two  consecutive  bits  of  RbI     /f        interval  between 

Generation  o-f  another  sequence  of  m 

are  in  practice  a  spreadTc    r^f      Y   nUmberE  whi<=h 

The  specific  sequence  used  i     ?'  S""d^  «*«»«. 

evetem  is  normally  not  avai  l-hf  I  com,T,uni"tions 

intended  receiver    (In  th  I    +  °  any0ne  but  the 
have  generated  the  sDreJn    Particular  simulation  we 
uniformly  distributed  noise'th^T*  by  P*5£ln9  * 
b't  interval  is  l//c  ,  9h  a  hard  "■"■!».   The 

Generation  of  thn   k 

?i-n  binary  .^^Tlffi^""  "PPln°  °<  th* 

-  -v.  into  either  .  or  .  ^^Vsa^ ?>« £. 
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DSSS  BPSK  Generator 


Uniform 
Random 
Number 


. 

RBS 

0.5 

X(i) 


Signal 

Cosine 

1 

X(i) 

Pul86 

r 

P- 

_Tk 

Expansion 

-^ 

> \> 

~T 

(1  •  kf  ) 

,. 

n 

2  f  k 

Kev 

1    0.5 

Figure  4-1.   Generation  of  the  BPSK  Signal 
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<See%TePtl°f   thD\fi:mt   hard  "-"t  is  stored  in 

the  array  y.   This  fnf^  ?*  the  Second  ^-d      e  ".f  f^  * 
4-1-         •    mis  information  1c  ^«+  •        1"llier  is  stored  in 

the  array  index  given  as  "IT  IL  *"   by  a  subtle  use  of 

discrete  sample  number  and  /  i"  th   J!  1S  the  lnde;;'  k  is  the 
intelligence.   Similarly  fnnfh     *  blt  rate  °*  the 

the  equations^*  ""lrt'»-   *■  ^ed^i  ^  J^ 


"%    -    Acos(2ir/.k    +    *    ) 


**    =   Cx(i)    0  y(j)2Tr 


c 

ow   is   given   by 

(4-5) 
(4-6) 


P  =  n 


Of  i-      = 


F 


k  /  y     k 


P=l 


(4-7) 


where   i    and    j    are   indices   n-f    fh 

^V'    i"P°rt*nt   point   to  note   thatThyS  **  defi?ed  earli^.      It 

^ined   by  the  parameter   ptiT//^1    u  *    ls  completely 

This   siqnal    is   not    rD.,i 
additional    feature  of   spreadLn   It""1*  BPSK  S19nal    but    it   ha=   an 
chip-frequency  and  carrTer   fr~n        B  Spect™»  by  control  lino   the 
diagram  of   the   scheme's  given  ^n"^  *""  baud"^te.      The        o 
V,    is  g1Ven   by:  "   9'Ven   ln   Fl9^  4-1.      The  desired  signal 


(4-8) 


Wher 


P-ts'of^heSfreguency°: pectr u"* V'9"*1    Pla"d  at   d'"-ent 
«gr«.I    generated^  I  ^Clflc  Xfk!"  ^  '^  b*"d  «* 


■J-3- 


4.4     Frequency  Shift   Keying   Sequence 

,h.c   -pouenre  needs   a  random  binary   intelligence 

distributed  noise  through  a  bar        »^^^t   i=k/fa   where 
f^the  raur-raayte'>o/the1n'.ormat1on   and    ,    is  discrete  sample 
number.    Now  the   desired   signal    is   generated   via: 


5|#    =   2cos(5k)sk_1    -   sk_2 


8k   =   B   +  £x(i) 


=  s„  +  r, 


(4-9) 

(4-10) 

(4-11) 


where  •  1.  the  -rier  freo^ ,  and  S  is  the  ^port'ant  and 
r^re'^nlucn-thl^r.tU  and  s.2=-sin<*,  giving  an  unit 

amplitude  sinusoidal  signal. 


4.5  Simulation 


chapters  will  now  be  s.«ul   «i  an  ^^ons!  These  simulations 
generated  as  described  in      f  q      °and  were  written  in 

retrCanr77ed  Fu  therVet  i   on Yhese  stations  can  be  found  Inl 
Fortran  ".   Furthe,  «  ^  are  the  Kwan  and 

^rfTfUter^tn  Modified  Kwan  anc <  Hart in  filter  new 
algorithm)  as  described  in  chapters  2  and  3.   The  ke>  P 
of "these  -filters  are: 

(a)  Sharpness  of  the  notch  filter  defined  by  pole       ' 

position  (r   -  1  ^       ,  .       nljlirB  /.,} 

(b)  Step  size  in  the  incrementation  procedure  (*) 

(c)  Time  constant  o-f  the  fading  filter  (X) 

(d)  Model  order  (n)  i  nfpr-f  erers  (m' 

(e)  Order  o-f  incoming  signal  or  number  o-f  interferes 
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4.5. 1      Fading   Filter 

In   order   to   tracJ    th      + ■ 
sssent.ally  achieved  by  a  simple^rajng ^ttr:    ^  ""  ** 

to  the  f,rst-order  forgrttfng^J^^}-  '-  equivalent 

-^rpVe^nu^^r  2^.^**  th'  *«*«  Cb,  and  Cc,  by 

°ptimi21ng  ,  and  X.   The'va  ue  o J*,  "M  ^^  by  ^nuaHy   " 

1UP  °f  *  was  cho^n  as  0.01  and  X=0.9. 

using'^rro^'Band^oi^'sTgn^l  v  ^  *  ■  ?"  aI9°ri^  was  tested 
signal  was  used  only  to  tune  th."  f"  ltS  P^^^ance.   This 
the  notch  filter,.   It  wL^  ^ha^if"  /J"'  ^pness^f 
clos»  unity  variance  of  the  naralf?    ,  i     the  Value  °*    r  is 

Based  on  simulation  the  value  of  f was  £  eStilMte  '"^ases. 

ue  °+  r  was  chosen  as  0.9. 

4-5.2  Check  VAX  verslon  agaanst  IBM-PC  version 

^theXtr::  ^iz^^rt^*^  *  -  •«.»- 

the  new  futer  that  win  be  £*£   ?r\?  °"S..the  block   dlagram  of 
notational  changes,  this  is  the  !L="   r*  ChaPter"   0ther  than 
shows  the  results  of  a  si ne-wave  tilt    ? ■  F ^^  ^   Figure  «-* 
chapter  3.   The  purpose  of  thU  tes J      !   ""'  descri°*d  in 
vers, on  of  the  program  gave  the  slm»     ^  Veri'y  that  the  VAX 
version  used  in  the  previous  «L!    reSUUs  as  the  IB«"PC 
insisting  of  three  sine  waves  fn  tT'  ^^   4"3a  was  the  input 

2  ♦  fK9RalS-  Fl^   4!3b  ls  t  e  fmrT^f  °f  n°1Se'  but  no 
that  the  filter  dad  remove  the  sfn!        °UtpUt'  de«onstr.ting 
parameter  adaptation.   By  ,   „  2   H  T^'   FlgUre  4"3c  sh°"*  the 
example  on  both  the  VAX  and  ne'lBM-PC   **  Comdaris°"  °*   this 
that  the  VAX  program  was  wort  i no  th/  '    HBre  able  to  verify 

working  the  same  as  the  IBM-PC  program 


■-'•_'~ 


Figure  4-2.   Block  Diagram  o-f  New  Algorithm. 
(Noil:   Same  as  Figure  2-4,  except  for  notation) 
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4.5.3     Elimination  of   narrow-band  BPSK  interference 

Our   next   experiment   was  ^signed  to  £  ^£  ™hm  sine 
algorithm  performed  or J^^^ftne  in  ut  consisting  of 
Wave   interferes       Figure  4  4a  s  ,requencies  as  the 

three  narrow-band  BPSK  "9na'=  *  shows  the  entered  output 

"^/^"that  Thernew"ftne        sable  to  eliminate  narrow-band 
Zr Tls  we  11    as  sine  lives.      Figure  4-4c   shows  the  parameter 
adaptation  for   this  case  and   should  be  compared  to  Flgure  4-3c. 

4.5.4     Transient  behavior    (tracking  a  moving  interfered 

The  transient  behavior   of   the  adaptive  filter  was  tested  by 
T     n  l  sianal    usina   the  hardware  of   Figure  4-1.      The 

generating  an   FSK  signal  the  5pectrum  ls 

FBK  signal    is  similar   to  the  Br»        8        ,  switched  back 

constantly  changing.      In  our   example,   ™e  SP  the  FFT  of 

and   forth   between   two  W-nci-^  F  9ure  ^sh   ^  ^^ 
the    input    signal.      From   hi  pure        j«  r:nlirp   4_«=;b 

fluencies  that   the   input   i^ang-g  ^«ee n.      F  gurej^b   ^ 
shows  the  output   specrU„  ^"^  the  parameter 

^irr^^Ln  t^alues       When -mpared  -^input 

^-rgThrFSr^lInd  Sl^CS  extremely  well    ««. 


4    5    ^ 


Elimination  of  narrow-band  BPSK  from  broad-band  BPSK 

Our  next  experiment  was  aimed  at  eliminating  ^ro>H>an^  | 

BPSK  interference  from  a  broad-band  BPSK  si gna 1.  figure  4-6a    . 
show-  the  input  consisting  of  three  narrow-band  BPSK  signals 
Bering  with  a  broad-band  BPSK  signa  J^P™^ 

with  the  desired  BPSK  signal,  we."e^,a"^°Jy",fecting  the 
aloorithm  eliminated  the  noise  without  adverse      usina  the 
signM  [1H.   Furthermore  -  P^^  ^  3  ^  'o "ained'the 

Cyclic  Spectral  Analysis  Software  Package  lioj 

Output  snown  in  Figure  ^7.   This  output  was  as  expctd  and    , 

confirmed  that  the  new  adaptive  notch  filters  were  p 

well  C123. 


-58- 


F:0ure4-4.   spectra,  of  3  Narrow-Band 


BPSK  signals  ln  Noise 
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Figure  4-5.   Tracking  of  an  FSK  Signal  in  Noise. 


-60- 


Input  . . 

Number  o 

Samples 

interfere 

Samples  p 

interfere 

Broad  ban 

carrier  c 

Broad  ban, 


chi 


P  bit 


Adaptive  filt_ 
Adaptive  filter 


f  interferer 
Der  carrier 
JS  .  . 

s  .  . 
cycle 

>er  Chip 
rs  .  . 

bit 

of 

d  signal 
ycle  . . . 

d  signal 

samples 

4 

samples 

er  step 


DS-BPSK 


of  the 
12.0,4.5,3. 

200,200,200 

per 
4.  0 

per 

.  0 
0.  02 
0.  9 


Figure  4-6.   Filtering  Narrow-Band  BPSK  «    D 

«**  Slgnal  in  whl t^S  ^™  B— ^nd 


-61- 


3.82M  -i 


a. 


Broad-Band  BPSK  without  Interference  or  Noise 


11. SMI  -1 


e  ee: 


b.   Broad-Band  BPSK  Corrupted  with  Interference  and  Noise 


3.821  -i 


c.   Corrupted  Output  After  Processing  by  New  Notch  Pilters 
Figure  4-7.   Processing  with  Cyclic  Spectral  Analysis  Software 
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4-5.6     Performance  with   too  many  notches 

One   situation   of    nrPaf    r-„„ 
there  were  too  many  notches  for  fT   "a\What  WOUld  hW»   if 
^gnals.      would  the  remainin „     t  mttber   °<   in^rfenng 
broad-band  BPSK  signaf?  - f   s   ^uld   1^^   tD  ^^  the 
Spectral    Analysis  Software  Package       n™  f  *"Ct   the  C^llc 

this  important    issue.      In   Figure  V-B  fh!     **  4"S  a"d  4"9  add^" 
and   three  notches  and   in  Figure  4-9  th^  tW°   inter'^ers 
and   three  notches.      As  before   V  *  ls  only  one  interferer 

spectrum,    Figure  4-8b  sn^ws  the  ^\4**  Sh°WS  the   '"<"* 
shows   the  parameter  adaptation        As'canT^1""    *"«  FlQure  4"8c 
the   remaining   notch   wanders   around   th,         *'*"    iR   R9ure   4-8c, 
lobe  of   the  broad-band  BPSK  signal    *h?/re7uenci*>  of   the  main 
However    Slnce   the   notch    is   very   Lrrn      P        9   t0   "nCel    it- 
of   the  broad-band  BPSK  siqna l      L2        ?"  ""Pared   to  the  main   lobe 
any   signif.cant   amount   o     s'onal    en"        "   '!  U°able  to  eliminate 
the  Cyclic  Spectral    Analyst  So!t!       ^'      Thu£  Wnen   we  applied 
results  as  Figure  4-7       b  a  =        , *  Packa9s,    we  got   the  same 

4-9-      The  two  remaining  notches TfxW^o  —  <°Und   in   F^e 
main    lobe  of   the  broad-band  BpIk  signa?     h   f    '   "^   ^^   'he 
eliminate  any  significant    am        1    signal,    but   are  unable  to 
Cvcl.c   Spectral    £.&   ^s™   e° ^al    energy.       Tnus   w°hen   ^ 
again   got   the  results  of   Figure  4-7   tlB?!  ^  3PPHed'    We  °nce 


4.6     Conclusions 
"he  BPSK  siqnal 


and  the  results   showed   ^hat'thTnew  iTl*"  ^  "ew  al^thm 
able  to  perform  well    under  various  sft^  "0tCh   fUters  were 
real    world.      The  new  aloori+h?  situations  expected   in   the 

mission  and  was  able  ^Zl^T'5^^  <*™«  «H    in   its 
°V  the  Cyclic   Spectral    An'afysTs'  S^™^*  <°'   P-cessing 
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Figure  4-8.   Two  Interferes  and  Three  Notches 
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5.0  Conclusions 


in  troad-band  communications  systems  in  the  presenu 
Gaussian  noise.   We  have  come  to  the  following  important 


conclusions: 
1. 


of   the  many  possible  algorithms  available  for  solving 
tnis  problem,  only  two  seem  to  meet  the  requirements 
ano  aLo  be  feasible  in  terms  of  hardware  complexity 
and  hardware  speed.   These  two  approaches  are:         j 

a    adaptive  second-order  notch  filters  based  on  a     I 
modified  Kwan-Hartin  approach  or 

b.   an  approach  based  on  the  DFT. 

This  report  has  studied  carefully  the  adaptive  notch   j 
'iter  approach  based  on  a  modification  of  hwan  and 
Martin  anu  found  it  to  be  extremely  well  suited  to  the 
problem. 

a.  Single  interfered  are  eliminated  within  at  most  a 
■few  hundred  iterations, 

b.  With  multiple  interferes  the  first  notch  will 
eliminate  one  of  the  interferes  wi thin  a  few 
hundred  iterations  and  other  notches  kick  in 
quickly  such  that  as  many  as  ten  interferes  can 
be  eliminated  in  as  few  as  50E  iterations,      f 

c.  When  there  are  fewer  notches  than  interferes, 
notches  can  be  designed  using  proper  *°r9e">"9 
filters  to  either  eliminate  the  maximum  number  of 
interferers  or  to  attempt  to  reduce  all 
interferers  while  eliminating  some  of  them, 

d.  When  there  are  more  notches  ^an  interferes,  the 
extra  notches  do  not  adversely  effect  the  broad 

band  signal , 
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e. 


i. 


togetherTand  °**  ^ferers  are  close 

hardware   impltL"TtToZ   ""  "" "  "eU    S^°   <°r 

5-1      Brief   Look   at    Hardware   Feasibility 

-  shall    look  br^lTat   pos si M.W   ^  ^  °<   ^  ^port, 
order   to  satisfy  ourselves  that  th=fG  conf  Rations   in 
desired  aDpllcation.        ^allv     J  !  *  9orithm   is  feasible  for   the 
notch   filter   that   woul  Toper Ite  I  , ' * "'  tD  deSi9n   *  5  to   10 

MHz.      With   currently  avanfbU  techrT        °   fre^uenci^  "P   to   100 
available  commercial! v  tnii*?        tffhnoloQy.    about   the  best 
Allows,    we  shall    a     eympt   to   Tuft^M"  '°  ™  MHz-      In   ** 
discuss  how  through   the  use  of  nfn   f-      "    10~20  HHz   fl?ure  and 
sampling   rate  to  as  much   al       g  ^  '"'"W1'   ^^  thl 
*ie:      80-160MHz).  9        tl(r'es  the  faasic   sampling  rate 

5-1.1      Commercially  Available  Hardware 

«"«.    .  shall    loo,:   ^STc^^r^^.^.^ 

''        aTl^l'e^idT^^!  CyC"  "d'    5  —  -«P1V. 

2-  3*s  sa  R-ri'^r^-irj;^41^ 


5 


-67- 


3.   Texas  .stents  TI  aB^  3J.  cycl.  2^1.^,  3 

cycle  multiply,  and  11  cycle  oiviu 

root)  C143, 

i    -j  rvrlp  add.  2  cycle  multiply, 
4    ATM  DSF32c,  20ns  cycle,  2  cycle  aaa ,  -  y 

and  3  cycle  divide  C153,  and 

5.    Honeywell  HDSP  6611B.  ^■^^i"^^^?? 

product  of  8  elements  -for  an  effective 

and  accumulate  time  C16D. 

for  the  time  being  on  the  other  four  chips. 

Figure  5-1  shows  the  blrtdlyj.*  ^^"^  The 

„„  th.n  b.  r,p..t.d  <=r  ,11  actlM  «~i  ™Jt.  J^^,  „, 
r.,»..r  .tnitur.  ••  '*'»«  'I^'^LJuJSi^lL.W  « 

calculating  the  filter  derivative  [103. 


5. 1.2  Time  Budget 

The  details  of  the  «»»>«*J^^^X££^ 

f^enTto  le^Lr  ESl'K/ElS  ^-—t.      However, 
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Figure  5-1.   Block  Dlagram  Qf  ^ 


d  Order  IIR  Filt 


er 
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Basic  Building  Block 
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Fig  5.3 


Figure   5-3.      Block   Di 


agram   for   3-Notch   Cancel 


er 
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this  was  done  without  pipelining  the  IIR  section   Using  the 
ninelinina  techniques  of  Soderstrand  and  Loomis  [17  21],  we  wouia 
h   sampling  rate  by  at  least  two  and  probably  four 
q  in  the  desired  10-20  MHz  design.   Furthermore  if  we 
were  to  vectorize  the  algorithm,  the  Honeywell  chip  could 
potentially  yield  100  MHz  sampling  with  both  vect°r^atl0?.^. 
^lining   Based  on  this  hardware  analysis,  we  feel  confident 

L,  rates  of  25  MHz  or  more  can  "-achieved  in  practice 
with  commercially  available  DSP  chips  or  most  certainly  wlth  a 
custom  VLSI  design. 


5.2  Future  Work 

There  are  two  avenues  of  future  work  that  should  be  pursued 
before  making  a  hardware  commitment  for  the  narrow-band  filter 

project: 

1    DFT  based  approaches  to  the  problem  should  be  studied, 
especially  in  view  of  the  availability  of  high-speed 
FFT  processors  and  such  exotic  approaches  at  recursive 
DFT's  using  QRNS  arithmetic  [93,  and 

a  detailed  investigation  o-f  hardware  implementations  of 
the  modified  Kwan-Martin  algorithm  should  be  made 
including  simulation  and  layout  of  a  custom  VLSI 
approach.   Because  of  the  regularity  of  the  structure 
tn  Figures  5-1  though  5-3,  the  custom  VLSI  may  well 
yield  some  very  excellent  filters. 
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